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The antiferromagnetic molecular wheel Feis of eighteen exchange- coupled Fe"'""'""'" ions has been 
studied by measurements of the magnetic torque, the magnetization, and the inelastic neutron 
scattering spectra. The combined data show that the low-temperature magnetism of Feis is very 
accurately described by the Neel-vector tunneling (NVT) scenario, as unfolded by semiclassical 
theory. In addition, the magnetic torque as a function of applied field exhibits oscillations that 
reflect the oscillations in the NVT tunnel splitting with field due to quantum phase interference. 



PACS numbers: 75.50.Xx, 33.15.Kr, 71.70.-d, 75.10.Jm 



Although magnetism is a priori quantum mechanical, 
observation of genuine quantum phenomena such as tun- 
neling and phase interference in magnets is difficult. This 
is because in ferro- or ferrimagnets, the low- temperature 
magnetism and dynamics are by and large well described 
by the magnetization vector M obeying classical equa- 
tions, although for magnetic molecules (e.g. Mni2, Feg), 
tunneling of the magnetization and phase interference 
have been observed [1]. Antiferromagnets, however, ex- 
hibit zero magnetization in the ground state, and obser- 
vation of quantum behavior is even more challenging. 

Antiferromagnets can be described by the Neel vec- 
tor n = (My^ — M^)/(2Mo), with sublattice magnetiza- 
tions M^, of length Mq (Fig. lb). In three dimen- 
sions, they exhibit long-range Neel order, but domains 
enable thermally activated or even quantum fluctuations 
of the Neel vector [2[. Nanosized, single-domain antifer- 
romagnetic (AFM) clusters would provide clean systems 
to study Neel-vector dynamics, but mono-dispersity is 
then a key issue. In small particles with weak magnetic 
anisotropy, the Neel vector can rotate [1, 1^ , while with 
strong anisotropy it is localized in distinguishable direc- 
tions, but fluctuates by quantum tunneling, i.e., Neel- 
vector tunneling (NVT) |J,@,0]- Typically, there are two 
tunneling paths (Fig. Ic) and interference occurs due to 
the phase of the wave function. This gives rise to char- 
acteristic oscillations in the tunnel splitting as function 
of applied magnetic fleld, which would be observable in 
static measurements such as magnetization [7]. 

Initial attempts to establish NVT with the ferritin 
protein [sl] unfortunately proved controversial owing to 
polydisperse samples and the presence of uncompen- 
sated magnetization Q . Recent attempts with the AFM 
molecular wheels CsFeg and Feio were encouraging steps 
forward [13, but were criticized with arguments that 



the tunneling actions Sq/H are too small and the NVT 
picture only approximately valid [uj. The latter at- 



tempts were stimulated by the theoretical prediction of 
NVT in such wheels [7] , and that its observation would be 
assisted by the monodispersity and crystallinity inherent 
in molecular compounds, and their result ingwell defined 
structural and magnetic parameters 0, [lH, [3, [3]. 

We here present measurements of the magnetic torque 
(r = M X B), magnetization, and inelastic neutron scat- 
tering (INS) on the AFM molecular wheel [Fei8(pd)i2- 
(pdH)i2(02CEt)6(N03)6](N03)6,orFei8 (Fig. la). The 
combined data show that the low-temperature mag- 
netism in Fei8 is very accurately described by the NVT 
scenario, as unfolded by semiclassical theory [7]. More- 
over, the torque is found to exhibit wiggles that directly 
reflect the oscillations in the tunnel splitting with fleld, 
or quantum phase interference indeed. 

The generic spin Hamiltonian for AFM wheels is 
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with nearest-neighbor Heisenberg interactions of strength 
J < (periodic boundaries assumed), a uniaxial single- 
ion magnetic anisotropy of strength D along the wheel 
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FIG. 1: (color online) (a) Molecular structure of Feis (red: 
Fe^^^). (b) Schematics of the classical ground-state configu- 
ration of the Fe spins (red, blue arrows) with the Neel vector 
indicated (green arrow, length not to scale), (c) Shape of the 
potential ^(n) in the high-field regime (q/jbD /huo = 1.2). 
The two tunneling paths from n = -\-y to —y are indicated. 



axis and a Zeeman interaction {N = 18 is the num- 
ber of ions, Si is the spin operator of the ith ion with 
spin s = 5/2 for Fe-"--^-^). The path-integral treatment - or 
semiclassical theory - of Eq. (1) has provided analytic re- 
sults for the low-T magnetism and describes its dynamics 
in terms of the intuitive picture of NVT . 

For D > and a field B along x perpendicular to 
the Neel- vector dynamics are governed by the poten- 
tial V(n) = N/\^J\{giiBBYnl ^ Ns'^Dnl, with hard- 
axis anisotropics due to the field \\x and the magnetic 
anisotropy \\z. Thus, V^(n) has minima at n = ±.y sepa- 
rated by an energy barrier (Fig. Ic). NVT between these 
two classical states occurs if the energy barrier is much 
larger than the attempt frequency hw^ = s^/8D\J\, pro- 
ducing a tunnel splitting A in the quantum-mechanical 
energy spectrum. In the high-field regime qiibB > hcoo 
the anisotropy barrier is smaller than that due to the 
field. Hence, the Neel vector tunnels via the z axis, 
through the anisotropy barrier U = Ns^D^ but the tun- 
neling from e.g. -\-y to —y may be either along the 'upper' 
or 'lower' path in Fig. Ic. However, the phase of the wave 
function acquires a term 7r/2 due to quantum fiuctua- 
tions plus a topological term iiN qiibB /\AJ\ proportional 
to the field. Thus, with increasing field one repeatedly 
tunes through destructive and constructive interference, 
and the tunnel splitting oscillates between zero and a 
maximum according to 
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with Aq = ShuJo^ySo/hex.^p{—So/h) and the tunneling 
action So/h = Ns^/2D/\J\. Sq/H is a crucial param- 
eter: the above condition for NVT to occur translates 
into Sq/H > 2 (henceforth we specifiy D/\J\ in terms 
of So/h). Importantly, the tunneling affects also the 
ground-state energy €o{B) = £{B) — A(5)/2, where 
£{B) varies smoothly, typically quadratically with B [?,]• 
Hence, the oscillations in A{B) can be detected, at zero 
temperature, via the static magnetization or torque [7]. 

Fei8 was synthesized as in Ref. [l^. Torque and mag- 
netization data were recorded on single crystals with ex- 
cellent quality. The shapes were three-sided pyramidal, 
allowing accurate alignment of the crystals. Torque was 
measured with a CuBe cantilever inserted into the MIO 
magnet at the Grenoble High Magnetic Field Laboratory 
equipped with an Oxford dilution fridge. Magnetic field 
was in the xz plane, with an angle 6 between field and z. 
Magnetization was measured with the same device, but 
displaced from the magnet's field center, which resulted 
in a Faraday force proportional to the magnetization. De- 
spite very careful alignment, a small torque contribution 
to the signal could not be avoided; hence only the parts 
of the data are shown where the estimated torque con- 
tribution is smaller than 5%. A smooth background is 
also present, which is estimated to be smaller than 5%. 
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FIG. 2: Magnetization vs field for (a) B\\z and (h) B ± 
and torque vs field for B nearly parallel to z as (c) measured 
and (d) simulated using Eq. (3) (T = 0.1 K). 



INS spectra of a polycrystalline sample of undeuterated 
Fei8 were recorded on the spectrometer IN5 at the In- 
stitut Laue-Langevin, Grenoble. Data were corrected for 
detector efficiency by a vanadium standard. Initial wave- 
length was A = 5.0 A and temperatures were T = 1.5 and 
4.2 K. Spectra were summed over all detector banks. Ex- 
perimental resolution at the elastic line was 95 /ieV. 

The magnetization and torque curves vs field measured 
at T = 0.1 K are shown in Figs. 2a and 2c for fields 
(nearly) parallel to z. Both curves exhibit staircase-like 
steps at regular field intervals, as is common in AFM 
wheels 0, [l2|, [l5|. However, for fields perpendicular to 

the magnetization, after a first broadened step, in- 
creases linearly with field, with only weak features ap- 
parent (Fig. 2b). The torque is even more striking: after 
a first step it exhibits wiggles at higher fields (Fig. 3a) - 
this is unprecedented for molecular magnetic clusters in 
general, and AFM wheels in particular [l6[. The wiggles 
in the torque directly correspond to oscillations in the 
NVT splitting {vide infra) ^ and their observation is the 
main result of this work. 

The sign of the torquejFigs. 2c, 3a) determines D > 0. 
As pointed out in Ref. [7|, the semiclassical theory then 
predicts a dichotomy in the magnetization as a key fea- 
ture of NVT, i.e., a conventional staircase-like behavior 
for parallel fields but a linear curve with weak features for 
perpendicular fields - exactly as observed in Feig. This is 
a first, strong indication of the accuracy of the NVT pic- 
ture in Fei8. The field derivative of the parallel magne- 
tization (Fig. 4a) and the positions of the maxima (steps 
in the magnetization) shown in Fig. 4b demonstrate the 
regular field intervals of the steps. The dependence of the 
first torque step on the field angle, Bi{0)^ is presented in 
Fig. 4c (features at higher fields show similar dependen- 
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FIG. 3: (color online) Torque vs field for B nearly perpen- 
dicular to z as (a) measured and (b) simulated {0 = 80.5°, 
T = 0.1 K). Panel (b) also shows the gap between the two 
lowest quantum states vs B as simulated with Eq. (3) (blue 
curve) and predicted by semiclassics, Eq. (2) (black curve). 

cies, but with much weaker variation). In AFM wheels 
studied before, the angle dependence is well described 
by Bi{6) = a + 6(cos^ 6> — 1/3), which was also observed 
for the wheels CsFeg and Feio with the largest So/h re- 
ported to date |18| . Fei8 shows a more pronounced 
angle dependence indicating a substantially larger So/h. 
Figure 4d presents the INS spectra at 1.5 and 4.2 K. 
Four transitions I, II, iii, and iv at ca. 0.3, 1.0, 0.8, and 
1.35 meV, and a weak feature v at ca. 0.6 meV are ob- 
served for positive energy transfer. At negative energy 
transfer, the expected anti-Stokes lines V and iv' are vis- 
ible. From the temperature dependence, peaks I and II 
are cold transitions, and iii, iv, and v are hot transitions. 

The combined data establish an AFM wheel, and 
Eq. (1) with a large Sq/H. The molecular structure of 
Fei8 nieans that not all exchange paths are identical, as 
assumed in Eq. (1). However, the energies of the low- 
lying states relevant in our experiments are not affected, 
to first order, by variations in the exchange constants 
fiot . J should thus be taken as the average, as is also the 
case for the anisotropy constant D. Thus, with a modula- 
tion in the exchange and anisotropy parameters, Eq. (1) 
is the appropriate model to describe the low-T behavior 
of Fei8 (and therewith also the models discussed next). 

Since the Hilbert space for Feis is huge 10^^), pre- 
venting direct use of Eq. (1), our quantitative analy- 
sis employs two approximations, which capture the low- 
energy i.e. low-T physics well. First, the AFM sub- 
lattices are replaced by spins = XIzga = 
X^iGB each of length Ns/2 = 45/2. This gives 

Hab = -JSa 'SB^d (^Sl, + Sl,^ + gfisS • B, (3) 



FIG. 4: (color online) (a) dM/dB vs field for B\\z (T = 0.1 K). 
(b) Fields of the maxima in dM/dB as measured (full circles) 
and simulated (open circles), (c) Field of the first torque step 
vs angle, Bi{0), as observed (circles), expected in the strong- 
exchange limit (dashed line), and simulated with Eq. (3) (red 
line), (d) INS spectra at 1.9 (full blue circles) and 4.2 K 
(open red circles). Blue solid and red dashed lines show the 
simulated spectra (elastic line not included in simulation). 

with j = aiJ and d = biD 0. The Hilbert space 
of Eq. (3) is small (2116), which permits exact numerical 
diagonalization. Second, we will use the analytical results 
of the semiclassical treatment of Eq. (1) [7]. 

Diagonalizing Eq. (3) and fitting to the data of 
Figs. 4b-4d gave j = -5.1(1) K,d = 0.021(1) K {g = 2.0). 
This parameter set reproduces all our data with high 
accuracy (± few percent), i.e., the fields of the mag- 
netization steps (Fig. 4b), the angle dependence Bi{0) 
(Fig. 4c), the INS spectra (Fig. 4d), and the torque curves 
for nearly parallel (Figs. 2c, 2d) and perpendicular fields 
(Fig. 3). Hence, Eq. (3), and therewith Eq. (1), describes 
the low-energy sector of Fei8 very accurately. The ob- 
served features in the torque (and magnetization) are 
less sharp than in the simulations (Fig. 3), and possible 
reasons for this include (i) additional, very weak terms in 
the spin Hamiltonian, e.g., Dzyaloshinski-Moriya interac- 
tions, which would lead to a rounding of the gap function 
and hence the torque, and (ii) a small distribution (5%) 
in the j and d values (j and d strain), which would result 
in a distribution of step positions and hence broadening. 

Determining J and D now requires knowing ai, bi. 
In principle, ai and hi should be determined such that 
the exact low-energy spectrum of Eqs. (1) and (3) match 
[lit . This is not trivial for Feig. Quantum Monte-Carlo 
techniques yielded ai = 0.2721 [i^l, and extrapolating bi 
for N = 6,8, 10 to Feis yields the estimate bi ^ 0.07. 
Hence, J = — 19 K and D ^ 0.3 K, which are typical 
values for ferric wheels. With these one gets So/h = 8.0. 
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We now turn to semiclassical theory. First a subtle 
point needs to be addressed. The above method of de- 
termining ai, bi is certainly correct, but we observe that 
semiclassical theory yields quantitatively accurate results 
only if one uses the ai, bi values predicted by itself. This 
is expected from self-consistency considerations. Hence, 
semiclassical formulae should be used with "corrected" 
J, D values as deduced with the semiclassical values 
af^ = 2/9, 6f = 0.1065 The tunneling action then 
becomes So/h= 5.90, which is the value to be used [22]. 

Henceforth we consider the high- field regime gjisB > 
hujQ, or B > 10.6 T for Feig. The criterion for NVT 
is Sq/H > 2, which is well fulfilled by So/h = 5.90. 
The expectation value of the Neel vector along ±y is 
{0|n^|0)2 ^ 1 - (5'o/^)~^ (3, hence in Feig the Neel 
vector is localized to 83%. Another figure of merit is 
the ratio Aq/U of maximal tunnel splitting vs barrier 
height. For Fei8, Aq = 0.320 K and /7 = 22.2 K, thus 
Aq/U = 0.014; the tunnel splitting is exponentially small 
as expected for a tunneling scenario. Figure 3b displays 
the energy gap between the two lowest levels for perpen- 
dicular fields as calculated with Eq. (3) and the semiclas- 
sical formula Eq. (2). The excellent quantitative agree- 
ment finally demonstrates the high accuracy of the semi- 
classical NVT theory for Feig. 

With reliable values for j, So/h, and the accuracy 
of our modeling of Fei8 established, we now discuss the 
wiggles in the torque for nearly perpendicular fields. Fig- 
ure 3b shows the calculated field dependencies for the 
energy gap and the torque. In the high-field regime, 
the energy gap (= tunnel splitting) exhibits the typical 
sinus-half-wave oscillation due to quantum interference 
0]. Comparison with the torque is striking: in the high- 
field regime the torque clearly follows the tunnel splitting, 
except for a smooth offset. The torque has not yet been 
calculated by semiclassical theory, but the result for Mx 
[7] inspires a refinement, which we checked numerically 
to describe the torque. For ~ 90°, we write the mag- 
netization as M = g jj^bB {Fx, 0, cos OFz) and the torque 
as r = g/j^sB cos 6{Fz — Fx), where, at K, 

Here we used Bq = 4\J\/ {NgfiB) and a factor (—1)"^ 
such that |sin()| = (— l)^sin() (n numbers magnetiza- 
tion/torque steps). As compared to the semiclassical 
formula for Mx we have, to grasp weak but important 
effects of magnetic anisotropy, introduced shifts <px{B), 
(j)z{B) to take into account that the steps do not occur 
exactly at fields nBo, and functions fx{B), fz{B) to ac- 
count for a smooth offset. Numerically we find that (j)x^ 
(j)z^ fxi and fz vary with field roughly as B~^ , in accord 
with expectation, such that after some rearrangements 



where g{B) and f{B) are smooth, essentially constant 
functions of B. The torque provides only indirect evi- 
dence for the size of the NVT splitting, but Eq. (5) proves 
that it does directly probe the oscillations in the NVT 
splitting, as demonstrated before in Fig. 3b. At non-zero 
temperatures, the sine term in Eq. (5) should be multi- 
plied by tanh(A/2/cBT), which has the effect of smearing 
out sharp features, or of rounding the torque wiggles. 
This is why the simulated torque curve in Fig. 3b is less 
sharp at the zeros than the gap. 

In summary, the combined magnetic torque, magneti- 
zation, and INS data, and their analyzes by two models, 
demonstrate that the low-temperature magnetism in Fei8 
is in accord with the NVT picture of semiclassical theory. 
Most importantly, the observed features in the magnetic 
torque have been demonstrated to directly refiect oscil- 
lations in the NVT splitting, i.e., quantum phase inter- 
ference. An interesting point is the number of electrons 
involved in the tunneling. In Fcig, the total spin on each 
AFM sub-lattice has a length of 45/2, i.e., 90 electrons 
are involved. Hence, the disconnectivity, which is a mea- 
sure for the " macroscopiness" of a quantum effect [23], is 
90. This makes Feig one of the most mesoscopic molec- 
ular systems exhibiting magnetic quantum tunneling. 
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